Which cosmological model with dark energy — 



phantom or ACDM? 



in 
o 
o 



Wlodzimierz Godlowski 
Astronomical Observatory, Jagiellonian University, Orla 111, 30-244 Krakow, Poland 

Marek Szydlowski 

Astronomical Observatory, Jagiellonian University, Orla 111, 30-244 Krakow, Poland and 
Complex Systems Research Centre, Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 

In cosmology many dramatically different scenarios with the past (big bang versus bounce) and in 
the future (de Sitter versus big rip) singularities are compatible with the present day observations. 
This difficulty is called the degeneracy problem. We use the Akaike and Bayesian information 
criteria of model selection to overcome this degeneracy and to determine a model with such a set of 
parameters which gives the most preferred fit to the SNIa data. We consider seven representative 
scenarios, namely: the CDM models with the cosmological constant, with topological defect, with 
phantom field, with bounce, with bouncing phantom field, with brane and model with the linear 
dynamical equation of state parameter. Applying the model selection information criteria we show 
that AIC indicates the flat phantom model while BIC indicates both flat phantom and flat ACDM 
models. Finally we conclude that the number of essential parameters chosen by dark energy models 
which are compared with SNIa data is two. 



I. INTRODUCTION 



The present observations of the distant supernovae type la indicate that the Universe is presently accelerating 0, Q 
due to the presence of some unknown form of energy violating the strong energy condition px + Spx > where px 
and px are energydensity and pressure of dark energy, respectively. While the different candidates for dark energy 
\ were pro p osed 1 ^ ,01 and confronted with observations |a, 0, 0> S IS El El , the cosmological constant A and phantom 



> 



■ fields |12|; Ilia LLiJ violating the weak energy condition px+Px > are most popular. Whereas both the cosmological 
constant and phantoms fields, described by the barotropic equation of state pi = wip (wi < —1), are negligible in 
the neighborhood of the initial singularity, they dominate late time evolution. In the case of the A domination we 
obtain the de Sitter state as a global attractor in the future and for phantom fields there are a big rip singularity. 
Both the singularities are generic features of their models. In the case of the big rip singularity at the some finite 
Q H i time the both scale factor and energy density are growing at infinity and the unexpected future singularity of density 
appears [TEL ITfj| . Alternatively the universe evolved from the Planck region where quantum gravity effects would 
have dominate. In this context quantum cosmology is used as a first approximation to the quantum gravity expected 
effects. Recently the proposal of quantum loop cosmology offers possibility of avoiding an initial singularity [l7llll|. 
For completeness we also consider the brane model ^| |2(j, |^| and the models with dynamical dark energy satis fying 
the dynamical equation of state px = wx(z)px, where wx(z) is linearized around z = (the present epoch) |2Sll23j. 
• i— i , Let us consider four different prototypes of dynamical behavior which take place in the neighborhood of the initial 
and final singularities. They give rise the seven representative scenarios, namely, the CDM model with A (ACDM), 
the CDM model with topological defect (TDCDM), the phantom CDM model (PhCDM), the bouncing ACDM model 
(BACDM), the bouncing phantom CDM model (BPhCDM), the brane ACDM model in Randall Sundrum version 
(Br ACDM), and the model with the dynamical equation of state parameter linearized around the present epoch 
(DEQS). They are completed in Table [I] where the dependence of the Hubble function on redshift z is given together 
with a number of models parameters. 

We consider the FRW dynamics in which dark energy is present. The acceleration of the Universe is due to 
the presence of dark energy for which the equation of state is px = wxp and wx = const < (quintessence) or 
px = wx(z)p and wx is variable with the cosmological time scale factor or redshift z. Therefore the basic equation 
determining the evolution is 

where p e s(a) is effective energy density of noninteracting "fluids", k = ±1,0 is the curvature index. Equation Q can 
be presented in terms of density parameters 
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TABLE I: The Hubble function versus redshift for the seven evolutional scenarios of the FRW models with dark energy. 
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where ^ = Q c fi{z) = O m ,o(l + z) 3 + Qx,of( z ) an d fim,o is the density parameter for the (baryonic and dark) 
matter scaling like a~ 3 . For a = ao (the present value of the scale factor) we obtain the constraint fi e ff,o + ^fc,o = 1- 
We assumed that energy density satisfies the conservation condition 

Pi = -3H(pi +pi), (3) 

for each component of the fluid p e g = T>pi. Then from eq. JSJ) we obtain the constraint relation q + fij^o = 1. In 
Table[l]we denote the numbers of parameters by d and also presented the names of the free parameters in the models. 
Please noted that for flat model, Slk,o = the number of the models parameters is equal d — 1. 



II. DISTANT SUPERNOVAE AS COSMOLOGICAL PROBES OF DARK ENERGY 



Distant type la supernovae surveys allowed us to find that the present Universe is accelerating |]J, |2j . Every year 
new SNIa enlarge available data sets with more distant objects and lower systematics errors. Riess et al. |24| compiled 
the latest samples which become the standard data sets of SNIa. One of them, the restricted "Gold" sample of 157 
SNIa, is used in our analysis. 

For the distant SNIa one can directly observe their the apparent magnitude m and redshift z. Because the absolute 
magnitude M of the supernovae is related to its absolute luminosity L, then the relation between the luminosity 
distance dz, and observed m and absolute magnitude M has the following form 

m-M = 51og 10 di+25 (4) 

Instead of df,, the dimensionless parameter Dl 

D l = H d L (5) 

is usually used and then eq. I@J changes to 

p = m- M = 5 log 10 D L + M (6) 

where 

M = -51og 10 i? + 25. (7) 

We know the absolute magnitude of SNIa from the light curve. The luminosity distance of supernova can be obtain 
as the function of redshift 

^^(H^^kfcf^) (8) 

where o — — 777 and 
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T{x) — sinh(x) for k < 

T{x) =x for k = (9) 

T{x) = sin(x) for k > 

Finally it is possible to probe dark energy which constitutes the main contribution to the matter content. It 
is assumed that supernovae measurements come with uncorrelated Gaussian errors and in this case the likelihood 
function C can be determined from the chi-square statistic C oc exp(— x 2 /2) where 

/..theor ( .obs^2 

X 2 =E ( ' \ l } . ( 10 ) 
i ° l 

while the probability density function of cosmological parameters is derived from Bayes' theorem pj. Therefore, we 
can perform the estimation of model parameters using the minimization procedure, based on the likelihood function. 
Especially we can also derive a one-dimensional probability distribution functions (PDFs) for some model parameters. 

For deeper statistical analysis we use the Akaike |25| and Bayesian information criteria |26|. It is crucial for our 
investigations to ask how well different theoretical models fits the same data set |27|. Because from the theoretical 
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TABLE II: The values of AIC and BIC for distinguished models (Table G} both for fiat and non-flat model. 
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FIG. 1: The one-dimensional probability density distributions (PDFs) for S7 m ,o, obtained for (from left to right) the ACDM, 
PhCDM, TDCDM models, respectively. 



investigation arising a large number of candidates for dark energy description it is necessary to select the essential 
model parameters from the considered set of parameters of the most general model. 
The Akaike information criterion (AIC) is defined in the following way pl| 

AIC = -2 ln£ + 2d (11) 

where C is the maximum likelihood and d is the number of the model parameters. The best model with a parameter 
set providing the preferred fit to the data is that minimizes the AIC. 

The Bayesian information criterion (BIC) introduced by Schwarz [2^ is defined as 

BIC = — 2 ln£ + din JV (12) 

where N is the number of data points used in the fit. 

The effectiveness of using this criteria in the current cosmological applications has been recently demonstrated 
by Liddle 27]. Taking CMB WMAP data [2^ and applying the information criteria, Liddle [13 found the number 
of essential cosmological parameters to be five. Moreover he obtained the important conclusion that spatially-flat 
models are statistically preferred to close models as it was indicated by the CMB WMAP analysis (their best-fit value 
is fitot.o = £A\o = L02 ± O- 02 at lcr level). 

Note that both information criteria have no absolute sense and only the relative value between differential models 
are physically interesting. For the BIC a difference of 2 is treated as a positive evidence (6 as a strong evidence) 
against the model with larger value of BIC 

The results of calculation AIC and BIC in the context of dark energy models are presented in Tableland Table ITTT1 
In Table HTl we show results for all models considerated for both flat and non flat cases without any assumed extra 
priors. In the general case the number of essential parameters in the cosmological models with dark energy is in 
principal two, i.e., (.ffoj^m.o)- It mean that flat model is favored in the light of the information criteria. We can 
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FIG. 2: The value of AIC in respect to fixed value f2 m ,o for three flat models (with topological defect, with cosmological 
constant and with phantom) with only one parameter Ho estimated. 



TABLE III: The values of AIC and BIC for distinguished models (Table Pi, with priors f2 m ,o = 0.3 both for flat and non-flat 
model. 
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FIG. 3: The value of BIC in respect to fixed value f2 m ,o for three flat models (with topological defect, with the cosmological 
constant and with phantom) with only one parameter Ho estimated. 



observe two rival models which minimizes the AIC and BIC. They are the ACDM model and phantom CDM (PhCDM) 
models. One can observe that both BIC and AIC values assume greater values for phantom models. This evidence 
can be regarded as a positive evidence in favor of the PhCDM model. 

We can derived one-dimensional probability distribution functions (PDFs) for the models parameters. In Fig. 
we present PDFs for Sl m o for all three flat models with two free parameters, i.e., for the flat ACDM, TDCDM and 
PhCDM models. There is a significant difference between predictions of these models. The ACDM model prefers a 
universe with fl m> o close to 0.3, the PhCDM model favors a high density universe while the TDCDM model favors a 
low density universe. In Fig.[3]and Fig.[j|]we present values of the AIC and BIC for these models. If Sl mj o < 0.22 then 
the information criteria favor the TDCDM model. For Q mfi G (0.22, 0.34), the ACDM is favored while for n m . > 0.34 
the PhCDM model is preferred. 

The similar analysis with the use of the information criteria is done in the case of the assumed prior Q m .o = 0.3 
31] (Table ITTTI) . In this case models with A is preferred over the models with phantoms, that is in contrary to the 
results obtained in the previous case of no priors for ^l m ,o- It clearly shows that more precise measurements of r2 mi o 
will give us the possibility to discriminate between the ACDM and PhCDM models. 
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III. CONCLUSION 



The main goal of this letter is to decide which model with dark energy is distinguish by statistical analysis of SNIa 
data. To do this we use the Akaike and Bayesian information criteria. The former criterion weighs in favor of the flat 
phantom model while the latter distinguishes the flat phantom and ACDM models. Assuming the prior Q m o = 0.3 
both AIC and BIC criteria weighs in favor of the flat ACDM model. 

It was considered the seven different models containing dark energy and dust matter (baryonic and dark matter). 
Our main result is that spatially flat phantom model is favored to ACDM model by information criteria of AIC and 
BIC. The further conclusions are the following. 

• The minimal number of essential parameters in the cosmological models with dark energy is in principal two, 
i.e., (i/oj^m.o)- The list of essential parameters may be longer, because some of them are not convincingly 
measured with present data, like the parameter w\. 

• The curvature density parameter docs not belong to the class of essential parameters when all the rest parameters 
are without any priors (with no fixed f2 m ,o)- At this point our result coincides with analogous result obtained 
by Liddle who found it basing on other observations. 

• If we consider models in which all model parameters are fitted then the PhCDM model with double initial and 
final singularities is distinguished. 

• When we consider the prior on f2 m ,o, then the model with two-dimensional topological defect has 51 m ,o < 0.22, 
the ACDM has Q m , G (0.22; 0.34), and the phantom model Q mfi > 0.34. 

To make the ultimate decision which model describes our Universe it is necessary to obtain the precise value of 
^m,o from independent observations. 

One can also draw from our analysis some kind of philosophical conclusion. Basing on simple and objective 
information criterion finally we obtain that SNIa data favored the models with the initial (big bang) and final (big 
rip) singularities. The big rip singularities are represents generic final future state of universe dominating by phantom 
dark energy. The bounce can be caused by high energy quantity modifications to FRW equations which make the 
cosmology nonsingular • If the energy density is so large then quantum gravity corrections are important at both 
the big bang and big rip. It is interesting that classical theory reveals itself boundness, i.e., classical singularities 
which can be distinguished by the information criteria. The account of quantum effects allows to avoid not only the 
initial singularity 0, but also escape from the future singularity [U |3 HH • 
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